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Abstract 

We report on a very accurate measurement of the static quark potential in SU(2) 
Yang-Mills theory in (2+1) dimensions in order to study the corrections to the linear 
behaviour. We perform numerical simulations at zero and finite temperature com- 
paring our results with the corrections given by the effective string picture in these 
two regimes. We also check for universal features discussing our results together with 
those recently published for the (2+l)-d Z(2) and SU(3) pure gauge theories. 



One of the most interesting open problems in lattice gauge theories is the construction 
of an effective string description of the static quark potential. Starting from the seminal 
papers of Liischer, Symanzik and Weisz [HE] much progress has been done in this direction. 
On the analytical side, the corrections to the linear term of the static quark potential 
induced by an effective string description have been calculated. These computations have 
been carried out considering effective string actions with the most general quartic self- 
interaction term and any type of boundary conditions for a rectangular L x R string 
world-sheet [3]. Recently a computation of the subleading correction has been performed 
also for the case of the three static quarks [I]. On the numerical side, the main prediction 
of the effective string theory, i.e. the presence in the static quark potential of the so called 
"Liischer term", has been tested in several pure gauge theories, ranging from the 3-d Z(2) 
gauge model to the 4-d S77(3) Yang-Mills theory [3J-[TT| . String effects have also been 
investigated in the excited states of the static quark potential: these studies indicate that 
an effective string description seems to hold when the distance between the quark and the 
antiquark is about 1.5-2/m [TOJ EH UH UHj • It is still to be understood how these results 
match with those of the ground state where the effective string description turns out to be 
valid already at about 0.5/m. 

In the last two years the development of new numerical methods - mainly the Liischer 
and Weisz's multilevel algorithm [18] - has allowed to test with much higher accuracy 
the effective string picture of the static quark potential. Furthermore, with the increasing 
precision of numerical results, new important issues concerning higher order terms in the 
effective string action can be addressed. 

Numerical studies have shown that at short distances and at high temperatures (but 
still in the confined phase) the naive picture of a free bosonic effective string is inadequate, 
however it is not clear how it should be modified. The simplest possible proposal, i.e a 
Nambu-Goto like action, seems to describe rather well the large distance/high temperature 
behaviour of the static quark potential [El HU [2H HH] , but it seems to be less successful in 
the short distance/low temperature regime ^31201 121] • Moreover, contrary to the Liischer 
term - which is universal and does not depend on the gauge group but only on the space- 
time dimensionality and on the string boundary conditions - there is no evidence (and no 
reason) that a similar universality holds also for the higher order corrections. It is thus 
interesting to understand to what extent these higher order terms depend on the gauge 
group of the theory. 

In this paper we investigate these questions by studying the SU(2) Yang-Mills theory 
in (2 + 1) dimensions. We have extracted the static quark potential from the 2-point cor- 
relation function of Polyakov loops both at very low (T < T c /7) and at high temperatures 
(T = 3T c /4). We have then compared our results in these two regimes with the predictions 
of the free bosonic string model and of the Nambu-Goto effective string model truncated 
at the second perturbative order. Moreover we have also compared our results with those 
obtained in [5] for the 577(3) Yang-Mills theory and in fUE] for the Z(2) gauge theory 



in (2+1) dimensions. In this study we have addressed the following three issues: 

• quantify the deviations from the Luscher term in the short distance/low temperature 
regime and test if they are compatible with a Nambu-Goto type action and/or to the 
presence of a boundary term. 

• study the large distance/high temperature regime of the static quark potential to 
check the reliability of the free bosonic string approximation as well as of the Nambu- 
Goto action. Notice that this is a non trivial test of the effective string description. 
In fact, in this regime, the effective string predictions can be obtained by a modular 
transformation of the short distance/low temperature results. 

• compare our results with the corrections measured in other pure gauge theories to 
study the presence of possible common features. 

This paper is organized as follows. In section Q] we discuss the main features of the 
effective string description of the static quark potential, listing a set of predictions which 
we shall later compare with the results of our numerical simulations. Section El is devoted 
to a general discussion of the (2+l)-d SU(2) Yang-Mills theory while in section El we give 
some details on the Luscher and Weisz's algorithm we have used to generate our data. 
Then, in the second part of the paper, we compare the theoretical predictions with our 
numerical data: first (section HJ) at short distance/low temperature and then (section |3J) in 
the large distance/high temperature regime. Finally, in section® we compare the Z(2), 
SU(2) and SU(3) results in the short distance/low temperature regime. Our conclusions 
are contained in section A preliminary account of the results presented here can be 
found in [161. 



1 Effective string predictions 

In this section we summarize some known results on the effective string description of the 
static quark potential. For a more detailed discussion we refer the interested reader to the 
original papers [TJ El or the recent review [TTJ . 

In a pure gauge theory the interaction between two static charges at temperature T = 
1/L (L being the lattice size in the periodic temporal direction) can be extracted by 
measuring the correlation function G(R) of two Polyakov loops P at distance R 

G(R) = (P(x)P(x + R)*) = e - F(R ' L) . (1) 



The free energy F(R, L) of the static quark-antiquark pair is then related to the static 
quark potential V(R) by 

V(R) = j-F{R,L) = -j-]ogG{R) (2) 

In the confined phase the force lines of the static quark-antiquark interaction are focused 
into a flux tube connecting the two charges. For not too short distances, the static potential 
rises linearly with a slope a given by the string tension: the leading behavior of F(R, L) 
is the so-called "area law" 

F(R,L)~aLR + k(L) , (3) 

where k(L) is a constant term depending only on L and with no physical meaning. In the 
strong coupling phase of the model - i.e. below the roughening transition - this behaviour 
can be obtained analytically by a strong coupling expansion. In the rough phase the 
transverse degrees of freedom of the flux tube become massless: the flux tube delocalizes 
and the above description represents only the leading behaviour. The suggestion of jT] is to 
try an effective description of the dynamics of the flux tube in terms of a fluctuating string. 
According to this effective string picture, in the rough phase of the theory one should modify 
cq.© in order to take into account the quantum fluctuations of the flux tube ("effective 
string corrections"). The exact form of these subleading corrections is unknown: they 
are expected to be a complicated function of R and L and they depend on the effective 
string action describing the flux tube dynamics. However, assuming that the string is 
smooth and the self-interactions are very weak and negligible in first approximation, one 
can consider the free bosonic string limit. This limit corresponds to take into account only 
the leading term of the perturbative derivative expansion of the effective string action, 
dropping out all the other terms depending on the flux tube dynamics. Hence the free 
bosonic string approximation amounts to consider a purely geometrical description of the 
flux tube. Indeed one finds that the corrections F^(R, L) to the area law behavior 

F(R, L) ~ oLR + k(L) + F?(R, L) (4) 

only depend on the shape and on the boundary conditions of the string world-sheet and on 
the number, (d — 2), of the transverse dimensions. In the particular case in which we are 
interested in this paper - i.e. periodic boundary conditions in the compactified temporal 
direction and fixed spatial boundary conditions along the two Polyakov loops - one finds: 

F?(R,L) = (d -2) log ( v (t)) ; -ir = ^ , (5) 

where r\ denotes the Dedekind eta function 

oo 

V(r) = q^l[(l-q n ) ; q = e 2 ™ , (6) 

n=l 



The labels q and 1 in F% recall this to be the first order term in the perturbative derivative 
expansion of the quantum fluctuations around the free string approximation. Eq. (J3J) is 
commonly referred to as the "free bosonic string approximation" . 

The Dedekind function has different expansions in the two regimes 2R < L and 2R > L, 
related to each other by the modular transformation r/(— 1/r) = yf—inqir). This yields to 
the following two different expressions for the quantum corrections 



2R<L 

2R> L 

F?(R,L) = 



+ ^log(l-e- 7mL/i? ) 



2AR 

n=l 



(d-2) , (7) 



" / '' + ^og^ + f>g(l-e- 4 ™*/ i ) 



6L 2 

71=1 



(d-2) 



The exponentially small corrections coming from the subleading terms in eq. (J7J) are neg- 
ligible unless we are in the intermediate region R ~ L/2. In eq. (JJJ) the leading term 
of the expansion is the well known "Luscher term" [T] which decreases with the distance 
as 1/R. Instead, in eq. (J8J) . the first term is proportional to R and to T = 1/L: it is a 
finite temperature correction which lowers the string tension as the temperature increases. 
Interestingly, the expression obtained for 2R > L, can also be reinterpreted as a low tem- 
perature result. If the temporal direction is looked at as a spatial direction and vice- versa, 
the correlator between the two Polyakov loops becomes the temporal evolution for a time R 
of a torelon of length L. The difference now is that the spatial boundary conditions of the 
fluctuating string are periodic and no longer fixed. With these different spatial boundary 
conditions the coefficient of the Luscher term turns out to be 4 times larger [3| as one can 
read out from the expansion (jHJ). 

In our study we want to compare the expectations coming from effective string descrip- 
tions of the static quark potential with the numerical results of Monte Carlo simulations. 
This comparison will be performed considering the two different regimes that we have 
discussed here above. In order to set the terminology, we refer to 2R < L as the short 
distance/low T regime and to 2R > L as the long distance/high T regime. Note that we 
are always in the confined phase of the pure gauge theory: hence, here, high T means that 
we are at a finite temperature not far below the temperature of the deconfinement phase 
transition. 

At high enough temperatures (i.e. for small values of L), higher order terms of the 
perturbative derivative expansion become important and they are no longer negligible. 
These terms encode the string self-interaction and depend on the particular choice of the 



effective string action. The simplest proposal - discussed for instance in (31 13 Ej _ is the 
Nambu-Goto action in which the string configurations are weighted proportionally to the 
world-sheet area. In this model the next-to-leading contribution to the free energy turns 
out to be (we have set d — 3) 

F q ,{R,L) = - Y ^^[2E 4 {r)-El{r)] . (9) 

where E 2 and E4 are the Eisenstein functions. They can be expressed in power series 
as follows 

00 

E 2 (r) = l-24j2*i(n)<l n ( 10 ) 

n=\ 
00 

E 4 {r) = 1 + 240 J2^(n)q n (11) 

n=l 

q = e 27TiT , (12) 

where ai{n) and cr 3 (n) are, respectively, the sum of all divisors of n and of their cubes (1 
and n are included in the sum). 

Besides the corrections coming from the self-interaction, there is also a boundary term 
related to the boundary conditions of the string. In fact, due to the presence of the Polyakov 
loops, the string has fixed ends and the effective string action may contain terms localized 
at the boundary. The simplest possible term of this type is [§] : 



A b = - dt 
4 Jo 



fdh\ 2 (dh 
\dz) z=0 \dz, ._ _„ 



(13) 



where b is a parameter with the dimensions of a length and h denotes the transverse 
displacement of the string; the factor | has been added to agree with the conventions of 
|Hj. This additional term can be treated in the framework of the zeta function regularization 
similarly to the computation of the free string correction. The result of this calculation 
(see (2111201 for details) is that, at first perturbative order in 6, the regularization of the free 
string action plus the boundary term gives the same result as the pure free string action 
(namely the Dedekind function discussed above), provided one replaces the distance R 
between the quarks by: 

R^R* = ?—- - (14) 

Thus, denoting by F^{R,L) the contribution to F q (R,L) due to the boundary term, we 
have at first order in b (remember that we have fixed d = 3) 

F*(R,L) +F*(R,L) = log r? fi^j (15) 



Exploiting again eq. © and the modular transformation of the Dedekind 77 function, this 
expression yields to the following behaviors 



2R* < L 



2R* > L 



F?(R, L) + F*{R, L) = -£L (l + A) (16) 



nR nb 1 , 2R b 

1 h - log 

6L 6L 2 5 L 2R 



F K R, L) + F*(R, L) = -^ + ^ + -Jo g ^-^ (17) 



Besides the static quark potential defined in eq. (J2J), the other observables we take into 
account are the force Q(R) 

Q(R) = V(R+1) - V(R) = -| log / G{n+ ' 



L & V G ( R ) 
and the combination |9 a c(R) related to the second derivative of V(R) 

c(R) = -±fr(V{R+l)+V{R-l)-2V{R)) = \fr±^ . (19) 

In order not to enhance the lattice artifacts when considering discretized derivatives of the 
static potential V(R), in these definitions we have considered the following quantities [221101 

R(R)- 1 = 2ir[A 2 (R-a)-A 2 (R)}/a (20) 

R(R)- 2 = 2tt [A 2 (R + a) + A 2 (R -a)- 2A 2 {R)} /a 2 (21) 

where A 2 (R) is the Green function between the origin and the point (R, 0) of the lattice 
Laplacian in 2 dimensions and a is the lattice spacing. In table Q] we report the values of 
R(R) and R(R) up to R = 20. Both quantities Q(R) and c(R) contain informations about 
the string effects in the static quark potential. Since Q(R) is related to the first derivative 
of V(R), it is easier to evaluate in numerical simulations but it depends explicitly on the 
string tension. Thus, if the uncertainty in the determination of a is not small enough, 
it can affect the subleading effects in which we are interested in. On the contrary c(R) 
only depends on the contributions related to the quantum fluctuations of the flux tube and 
can directly probe the reliability of the effective string description. If the effective string 
corrections are absent (like in the strong coupling phase below the roughening transition) 
then c(R) = 0. Hence, even if the numerical estimates for c(R) are in general less precise 
than those for Q(R), in the following we shall mainly use c(R) in our analysis. The only 
exception will be the study of the long distance/high T data in section where we shall be 
able to extract important informations also from Q(R) by using the high precision estimate 
of a obtained at low temperature. 
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4.432 
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5.448 
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6.458 


6.920 
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7.464 


7.932 
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8.469 


8.941 


10 


9.473 


9.948 


11 


10.475 


10.953 


12 


11.478 


11.957 


13 


12.480 


12.961 


14 


13.481 


13.964 


15 


14.483 


14.966 


16 


15.484 


15.968 


17 


16.485 


16.970 


18 


17.486 


17.972 


19 


18.486 


18.973 


20 


19.487 


19.975 



Table 1: Values of R(R) and R(R) (see eq.s J2DJ) and fSJ) up to R = 20. 



It is useful to write explicitly the predictions of the effective string description in the 
short and large distance limits for Q(R) and c(R). 

Short distance/low T regime: 2R < L. 

In the free bosonic string case, when 2R < L the Dedekind function can be approximated to 
—7iL/24R (i.e the Liischer term only), neglecting the exponentially decreasing corrections. 
In this way, we obtain 



Qi(R) 



71 
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and 



ci(R) 



R*_1T_ 

2 24 ( 



2VR+1 



2 1 

R + R+l 



B/ ~ a+ 2AR(R + 1) 
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Ti 



R? 



(22) 



(23) 



R-y 2AR(R 2 -1) 

Hence the free bosonic string action predicts that, as R goes large (always fulfilling the 
constraint 2R < L) , c(R) approaches the asymptotic value lim^oo c(R) = |j which is the 
well known coefficient of the Liischer term in 3-d. 



Let us now look at the corrections given by the Nambu-Goto string model truncated 
at the second order. For 2R < L, the Eisenstein functions can be approximated to 1 (see 



eq.s (fTUj) and (fTTj) ) up to exponentially small corrections and we obtain 



q^)=q^)~ ( ~trtw-^ (24) 



Similarly one finds: 



e 2 (fl)= Cl ( fl) + ^^(-A + _l_ + _l_ ) (25) 

Finally, if a boundary term is present, we find an additional contribution to be added to 

the previous ones 

bn 2R+1 

QbiR) = 24 R 2 (R+l) 2 ^ 

and _ 

c >< fl > = 12 (fl^Tp < 27 > 

Large distance/high T regime: 2i? > L. 

Similarly to the short distance/low T regime, we consider the free string approximation 
and the Nambu-Goto model. In order to construct the expansions in the limit 2R > L, 
one performs a modular transformation of the Dedekind and Eisenstein functions. For the 
Dedekind function this transformation is given in eq.(JHJ), for the Eisenstein functions they 
have the following expressions 

^..L. 4R2 ,.2R. 12R AR 2 

E 2 (i — ) = 7rE 2 {i — ) H 5- (28) 

2V 2R J L 2 2V L J irL L 2 y ' 

^..L. 16i2* ,.2R. 12R 16i? 4 

Neglecting again exponentially decreasing terms, we obtain the following expressions (as 
above, the indices '1' and '2' refer to the free string approximation and to the Nambu-Goto 
model respectively) 

Q l{R) = a- ^ + ± log(l + i) (30) 

ft(B) = CMR) - 4s ~ ^Rdht, < 31 > 

* (iJ) = Tlr(WTT) (32) 

and _ 

ci(i?) = -^log(l--^) (33) 



These are the expressions with which we shall compare our results in section [5J 

It is useful to expand the previous expression in 1/R in order to better understand the 
meaning of these terms 

„, n s ( 7T 7T 2 \ 1 1/6 1 1 \ 

«( fl »~^-6^-72ljJ + 2lS + ip(2l-8^-ilJ + ''' (36) 

and 

C W~5 + (^-8^) + "' < 37 » 

Eq. ()36|) shows that the effective string description can also account for finite temperature 
corrections to the zero temperature string tension. Eq. (|37j) instead predicts a linear 
dependence of c(R) with the distance R. 



2 The model and the simulation parameters 

We have studied the (2+1) dimensional SU(2) Yang-Mills theory on the lattice with the 
Wilson action 

S = PY,(l-h±Ua) (38) 

□ 

The sum is over all the plaquettes of a cubic lattice with L spacings in the temporal 
direction and N s in the two spatial ones. The gauge coupling is denoted by (3 and U a is the 
path ordered product of the gauge field along the plaquette. We impose periodic boundary 
conditions in order to consider the Yang-Mills theory also at finite temperature. 

The Yang-Mills theory in (2+1) dimensions is not scale invariant at the classical level 
since the continuum gauge coupling g 2 is a dimensionful quantity. It has the dimensions 
of a mass and it is related to the dimensionless gauge coupling f3 by 

\ = a/3 (39) 

where a is the lattice spacing. Hence, in (2+1) dimensions, scale invariance is explicitly 
broken in the Lagrangian and it does not result - as in (3+1) dimensions - from the 
dynamical phenomenon of dimensional transmutation. 



Close to the continuum limit all dimensionful quantities like the string tension at zero 
temperature <r(0) or the deconfinement temperature T c can be written as power series in 
ag 2 /4 = (3. The first few coefficients of these series are known quite accurately E 



^ = 1.324(12) + 1J0(11) 4 



aT a 



P 
1.50(2) 



3 2 



+ 



(40) 
(41) 



It is important to stress that the string tension u{T) depends on the temperature: it 
decreases as the temperature increases and vanishes when T — > T c . There are by now 
clear numerical evidences that the deconfinement transition is second order. Then the 
Svetitsky-Yaffe conjecture [21j suggests that the universality class of this transition is the 
same as that of the 2 dimensional Ising model. Thus the expected critical behaviour of 
the string tension is cr(T) oc (1 — T/T c ) 2u with v — 1. Many numerical investigations have 
confirmed with high accuracy the validity of this expectation. In the following we shall 
mainly use the zero temperature value of the string tension and we shall simply denote it 
by a. 

Our numerical simulations have been carried out mostly at j3 = 9, corresponding to the 
lattice spacing a ~ 0.074/m. At this value of the gauge coupling the inverse deconfinement 
temperature 1/T C is about 6 in lattice units and the thickness of the flux tube connecting 
the quark-quark pair roughly corresponds to R c = a/1.5/ct ~ 7.5. Hence R c can be 
considered as the scale below which the effective string picture is no longer expected to be 
valid. In order to investigate the scaling behaviour, we have also performed a numerical 
simulation at (3 = 7.5, corresponding to the lattice spacing a ~ 0.09/m. In tabled we have 
summarized the simulation parameters considered in our study. 



p 


L 


N s 


9.0 


8 


120 


9.0 


42 


42 


9.0 


48,54,60 


36 


7.5 


48 


32 



Table 2: Set of couplings (3 and lattice sizes N 2 x L we have considered in our study. 



We have studied the model in two distinct regimes: 



First we have studied the short distance/low T regime, considering lattices with a 
large temporal extent: we have chosen (3 = 9, L 



42,48,54 and 60. In order to 
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keep small the finite temperature string corrections - depending on the ratio R/L 
- we have measured the static quark potential up to distances R = 13. The lattice 
size N s in the spatial directions has then been fixed so as to make negligible the echo 
effects due to the periodic boundary conditions. We have set N s = 42 for L = 42 
and N s = 36 for the other three cases L = 48, 54, 60. The chosen set of parameters 
leads to several simplifications in the effective string predictions, allowing to measure 
accurately the Luscher term and to detect possible higher order corrections. Having 
results at four different values of L, we have also been able to study the dependence 
on L of the various corrections. In this regime we have performed a simulation at 
j3 = 7.5 on a 32 2 x 48 lattice to investigate the corrections to scaling. 

Second we have investigated the long distance/high T regime: we have considered 
L = 8, corresponding to a temperature of 3/4T c . In this case we have measured the 
static quark potential at much larger distances - up to R = 20 - in order to reach 
values of the ratio 2R/L much larger than one. Since the temperature is high, the 
value of the string tension is quite small. We have carried out our simulations on 
a lattice with spatial extent N s = 120 in order to make the echo effects due to the 
periodic boundary conditions negligible at R = 20. This choice of the lattice size 
has allowed us to explore the "modular transformed" region of the effective string 
predictions. In this way a set of very stringent tests both on the leading correction 
and on the possible higher order terms have been possible. 



3 The algorithm 

In this section we describe the numerical techniques we have used in our Monte Carlo 
simulations. Since the effective string picture is a description holding in the low-energy 
regime, one has to measure the correlation function of pairs of Polyakov loops very far 
apart. This is a hard numerical task due to the exponential decrease of the correlation. 
The algorithm recently proposed by Luscher and Weisz jTH] has been a breakthrough in the 
numerical estimation of correlation functions of Polyakov loops: it yields to an exponential 
reduction of the error bars w.r.t. the numerical methods previously used. 



The idea of the algorithm - in some respects similar to that of the multihit [25] - is to 
reduce the short wavelength fluctuations. By freezing the Monte Carlo dynamics of a suit- 
ably chosen set of lattice links, one splits the lattice into sublattices non communicating 
among themselves. Then the observables are built up combining independent measure- 
ments performed within every sublattice. Suppose that a measurement of an observable 
O is obtained using the results of averages O su b computed in Af different sublattices. If iV 
sublattice measurements have been carried out, the combination of the sublattice averages 
O su b corresponds to (N)^ measurements of O. Hence, at the cost of about N measure- 
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ments in the whole lattice, one gets an estimate of O as if (N)^ measurements would 
have been carried out. However, this estimate is biased by the background field due to the 
frozen links. The next step is then to average over the background field, computing biased 
estimates for many values of the background field. 



Let us consider the correlation function of two Polyakov loops 



(P(0)P(xY) = - hjdUy^TT 



(0,0),4 ' • • U (0,L-1),4 



Tr 



* * 

U{x,Q),i ■ ■ ■ U{x,L-l),4 



-S[U] 



U-l' 



(42) 



We now split the lattice along the temporal direction into M = L/rit sublattices with a 

(s) 

temporal thickness of rit lattice spacings. If we now keep fixed the set V fc of all the spatial 
links with time coordinates kn t , k = 0, . . . , (Af — 1), then the dynamics of every sublattice 
depends only on the background field of the two frozen time-slices that sandwich it along 
the temporal direction. Hence the sublattices are totally independent among themselves. 
Keeping this slicing in mind, we rewrite eq. (|42j) as follows 



(P(O)P(f )*> = f J] <H7« T^[V^, v} a) ] . . . T^[V^_ 1} , V (s) ] V\vt s) ] 

J 7„ 



(43) 



r(«) tA») 



M 



r( S ) 



where T[V^ , V ] is the following sublattice average with fixed spatial links V { and 



V- at the two boundaries in the temporal direction 



rrpa7/3<5 r-rr (s) T/0)l 
(d,x)L V i > \j J 



n dU M 



(0,0),4 • " " U (0,n t -1),4 



U4' 



a-y L 



"(x,0)A • • • ^(x,nt-l),4 



e -S[U;V[ B) ,VJ s) ] 

PS Z[v} s \vf s) ] 

(44) 



by 



The sublattice partition function with fixed temporal boundaries Z[V^ , V ] is defined 



W T/( s )l 






(fl/^e-^^^ 



'y.M 1 



(45) 



where S[U; V { , V^ ] is the Wilson action of the sublattice with fixed temporal bound- 

(s) (s) 

aries VJ and Vy . In eq. Q43j) a, (3,^,5 are color indices and the contraction rule is 
qpa7/3<5'j7(T5T _ jacr/jT^ j^. j g important to note that T a7/M are gauge-invariant quantities 
under sublattice gauge transformations. Finally the quantity V\Yl ] is the probability for 



the spatial links with time coordinates kn t , k = 0, . . . , (A/" — 1), to be V h 



(s) 
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nvt s) ] = UU dU y „ II ^ - V t ] ) e- S[U] (46) 

y,H k 

Eq. P^jl states that if we evaluate T Q7/35 in every sublattice by performing N sublattice 
measurements, then we obtain an estimate of (P(0)P(x)*) - at fixed V^ - as if (N) 
measurements would have been carried out. Then, changing the background field according 
to the probability distribution of eq. (pltjj) and repeating the above procedure, one completes 
the numerical computation of the integral of eq. (|43|) . The described technique is the 
so-called "single level algorithm". Luscher and Weisz have also presented a generalized 

(s) 

"multilevel algorithm" in which the updating frequency of the background field V fc is not 
the same for the various time slices k. However, it seems that the "single level algorithm" 
is more efficient [SJ and this is also what turns out from our experience. 

When performing a numerical simulation with the single level algorithm, one has to fix 
the value of 3 parameters: the temporal sublattice thickness n t , the number N of sublattice 
measurements of T/^-s and the number M of background field configurations to carry out 
the integration over V^ s \ These 3 parameters are related among themselves and finding 
their optimal choice is a not trivial step. Moreover, they also depend on the temporal 
extension L of the whole lattice and on the distance R = \x\ between the two Polyakov 
loops. Some results on the optimization step can be found in [T2J |2H] • 

The temporal sublattice thickness n t has a minimal value coming from the knowledge 
we have that the Polyakov loop correlation function decreases exponentially (P(O)P(x)*) ~ 
e~ aLR . If the sublattice temporal extension nt is such that the sublattice is in the confined 
phase, then TV^ ~ e"™ tfi . Hence we evaluate the small number e~ aLR as the product of 
J\f much larger numbers e~ antR . In this respect, the Luscher and Weisz algorithm exploits 
a similar solution the snake algorithm [27JI2B] uses to compute accurately the exponentially 
small value of the 't Hooft loop. An estimate of the sublattice minimal thickness is L crit /2, 
where L crit is the temporal extension at which the deconfinement transition takes place at 
a given gauge coupling j3. 

The final error bar of the measure of (P(0)P(x)*) is the combination of the uncertainties 
of the sublattice averages ^@g\ (depending on N) and of the fluctuations of T,^ due to 
different background fields (depending on M). If we suppose to have fixed n t , the larger is 
R the bigger must be both N and M. Typically, iV is order of several thousands and M of 
few hundreds. We note that N does not depend on L while M does since it is related to 
the number A^ of frozen time-slices. According to our experience, when 2L/L cr i t is large 
- i.e. at very low temperature - it is often better to consider sublattices with a thickness 
n t bigger than the minimal one. Although iV must then be increased, the reduction of the 
number of frozen time-slices M - and, hence, also of M - makes the choice convenient. 
Indeed, in our short distance/low T simulations at (3 = 9, we have considered n t = 6 even 
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if the minimal thickness was n* = 3. However, we have not performed a systematic study 
on the tuning of the parameters. 

As a last remark, we note that the Luscher and Weisz algorithm allows an exponen- 
tial gain in the accuracy of the numerical estimation of (P(O)P(x)*) only in the temporal 
direction. Indeed, every sublattice estimate T^) ~ e~ antR is a quantity decreasing expo- 
nentially fast with R but it is still estimated by "brute force", i.e. with an error reduction 
proportional to 1/x/iV. It is possible to devise a further, partial gain also in the measure- 
ment of the sublattice averages. Similarly to the above described time slicing, one can cut 
every sublattice by freezing space-slices so that every spatial direction is halved. Hence 
a d + 1 dimensional sublattice will be split in 2 d independent sub-sublattices. Then, if 
the spatial positions, and x, of the two Polyakov loops are in different sub-sublattices, 
it follows that K measurements performed in each single sub-sublattice will result in K 2 
measurements of the correlator. In the numerical results we present in this paper, we have 
not implemented this procedure for an additional partial gain: however it can become 
convenient when the two Polyakov loops are very far apart. 



4 Analysis of the short distance/low T data 

In this section we present the results we have obtained from the numerical simulations 
in the short distance/low T regime. The values that we obtained for Q(R) and c(R) are 
reported in tables El and |U respectively. 



R 


L = 42 


L-48 


L = 54 


L = 60 


2 


0.040462(11) 


0.0404558(87) 


0.040437(10) 


0.0404552(73) 


3 


0.034381(14) 


0.034373(11) 


0.034349(13) 


0.034374(10) 


4 


0.031433(17) 


0.031422(14) 


0.031394(16) 


0.031426(12) 


5 


0.029790(20) 


0.029778(16) 


0.029747(19) 


0.029785(14) 


6 


0.028779(22) 


0.028766(18) 


0.028733(21) 


0.028777(16) 


7 


0.028112(25) 


0.028099(20) 


0.028063(24) 


0.028113(17) 


8 


0.027648(27) 


0.027636(22) 


0.027597(26) 


0.027654(19) 


9 


0.027312(29) 


0.027302(23) 


0.027260(29) 


0.027322(22) 


10 


0.027061(31) 


0.027052(25) 


0.027009(32) 


0.027076(24) 


11 


0.026869(34) 


0.026862(28) 


0.026819(35) 


0.026887(28) 


12 


0.026717(36) 


0.026714(32) 


0.026672(40) 


0.026739(33) 


13 


0.026597(40) 


0.026600(38) 


0.026554(48) 


0.026616(40) 



Table 3: Values of Q(R) for various values of L at (3 — 9. 
It is interesting to compare our data with those of jT5j which were obtained in a 
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R 


L = 42 


L = 48 


L = 54 


L = 60 


3 


0.067293(40) 


0.067319(32) 


0.067382(39) 


0.067299(28) 


4 


0.083360(93) 


0.083418(74) 


0.083527(91) 


0.083345(66) 


5 


0.09519(18) 


0.09527(14) 


0.09543(17) 


0.09508(13) 


6 


0.10392(30) 


0.10402(24) 


0.10427(30) 


0.10367(24) 


7 


0.11058(48) 


0.11062(39) 


0.11104(50) 


0.10996(39) 


8 


0.11583(73) 


0.11558(63) 


0.11629(81) 


0.11465(64) 


9 


0.1198(11) 


0.1195(10) 


0.1205(13) 


0.1185(11) 


10 


0.1236(16) 


0.1226(16) 


0.1233(21) 


0.1212(18) 


11 


0.1266(23) 


0.1252(28) 


0.1250(35) 


0.1244(31) 


12 


0.1298(34) 


0.1262(51) 


0.1256(62) 


0.1262(56) 


13 


0.1310(54) 


0.124(10) 


0.128(12) 


0.134(11) 



Table 4: Values of c(R) for various values of L at f3 = 9. 

similar range of temperature and interquark distances. In particular, the sample which 
allows the most natural comparison is the one at /3 — 10 and L = 48 reported in table 2 
of J2]- While Q(R) contains a residual scale dependence due to cr, the function c(R) can 
be directly compared. It is easy to see that our values are fully compatible with the ones 
of ^2J and slightly more precise. 

The data for c(R) collected in table E] clearly show a smooth convergence toward the 



expected free bosonic string value, 
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0.13089...: the last two values are already compati- 



ble with this result. Furthermore, it also turns out that the data corresponding to different 
extensions L in the temporal direction are compatible. This indicates that - within the 
errors bars - not only the Liischer term, but also the higher order corrections to the free 
energy are proportional to L (so that they become L independent in Q(R) and c(R)). 

Based on these results, we have explored the higher order corrections studying the 
deviations of the numerical data of c(R) from the free bosonic string predictions of eq. (|5|). 
Thus we have considered the following difference 



A(R) ee c(R) + ^ (F?(R + 1, L) + F*(R - 1, L) - 2F?(R, Q) 



(47) 



The data for the four values of L are plotted in figure H As anticipated above, it clearly 
turns out that the data for different L nicely agree among themselves. It is important to 
address the question of the range of interquark distances in which we must expect these 
higher order corrections. As mentioned in section |21 there is a natural scale jllj . R c , below 
which one cannot trust any kind of effective string description. In fact at this scale the 
interquark distance becomes comparable with the string width: the internal structure of 
the flux tube becomes relevant and the approximation with a string is no longer valid. 
Below this scale, one actually enters in the perturbative domain of the theory (see the 
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discussion in and [121 )• At the value of the gauge coupling we have used in our numerical 
simulations, we have R c ~ 7.5 and so we must consider only distances R > 8. Moreover, 
as discussed above, the numerical data at R = 12, 13 are already compatible with the free 
bosonic string approximation and so the higher order corrections are zero within the errors. 
Thus we concentrate our analysis on the four values R — 8,9, 10, 11. Our first goal is to 
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Figure 1: Deviation of the results of our numerical simulations at L = 42, 48,54, 60 and 
of the Nambu-Goto string model truncated at the second order (continuous line) from the 
free bosonic string approximation. 



estimate the dependence on R of the correction to the free bosonic string behaviour. To 
this end the simplest option would be to fit the data. However with this approach one 
has to face a major objection. Due to the numerical method we use, data at different 
distances R are correlated and the crosscorrelation matrix is so flat that keeping it into 
account in the fits is very difficult. To circumvent these problems we decided to assume 
a given correction, say b a /R a , extract the coefficient b a from each one of the four data at 
our disposal and then check the stability of the results obtained in this way. Let us see 
this in more detail. Let us assume that the potential is described by 



V(R) 



j-G{R) 



„ , F?(R,L) b a 
L R a 



(48) 
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then it is easy to see that A(R) should behave as 



A(R) 






+ 



(R + l) a (R-l 



2 
R~ a 



(49) 



From this equation one can estimate the value of the coefficient b a . It is interesting to see 
the physical meaning of this parameter. A typical string self-interaction term is modeled 
by the insertion in the Lagrangian of a quartic interaction term (as in the case of the 
Nambu-Goto action) and leads (see for instance [3] or ^T]) to a correction with a = 3. 
Instead a boundary term would lead to a correction with a = 2. In particular, with our 
choice of notations, we have b 2 = — ff (see eq. lfTfil) ). 

As an example of the results of our analysis we report in table 03 the values of b a for 
a = 2 and 3 for the data collected at L = 42. For the case a = 3 the estimates at 
the various distances R turn out to agree within the errors. The results obtained using 
the data collected with the other values of L are fully compatible with those reported in 
table Combining together the results collected at the four different values of L (they are 
completely uncorrelated since they have been obtained from independent simulations), we 
obtain as our final estimate for 6 3 the value 63 = 0.145(15) 



R 


b 2 


h 


8 
9 
10 
11 


0.038(2) 
0.032(3) 
0.022(5) 
0.014(9) 


0.15(1) 
0.14(1) 
0.11(3) 
0.08(5) 



Table 5: Values ofb 2 and b 3 from eq. XW\) using the data collected with L = 42. 

Since the parameters b a are dimensionful quantities, it is useful to study the dimen- 
sionless ratios 

la = b a a<>*-W (50) 

This will allow us to compare our results with other pure gauge theories, with other SU(2) 
simulations performed at different couplings (3 or with a theoretical model like the Nambu- 
Goto one. For instance, in the a = 3 case, we find: 



73 = 0.0038(4) 



(51) 



Unfortunately, the data do not allow to disentangle unambiguously a boundary term 
contribution from the one due to a quartic correction in the string action. However the 
stability of the numerical results for b 3 suggests that - similarly to the 577(3) case in (2 + 1) 
dimensions - there is no boundary term. Nevertheless, if we try to interpret the deviations 
from the free bosonic string as a boundary correction, we find that the boundary parameter 
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b must be negative and rather small, ranging from b ~ —0.3 for R = 8 to b ~ —0.1 for 
R — 11. Further investigations are needed to give more definite statements in this respect. 

Instead, if we assume that the boundary term is absent and that the correction is 
completely due to a selfinteraction term in the effective string action, then from our analysis 
we obtain two important results which we present here and shall discuss in greater detail 
in section EJ 

1] The value of 73 that we find is similar (but not equal within the errors) to the one 
obtained in the Z(2) gauge theory at (3 = 0.73107: j 3jZ(2) = 0.0047(1) 



2] The value of 73 that we obtain is definitely different from the Nambu-Goto expectation 
7 37VG = — YY52 — —0.00857.. (see eq.® ) which turns out to be opposite in sign and 
more or less double in magnitude. It is important to stress that a similar disagreement 
in the short distance regime was recently reported in ^UJ E3 |!H] . 



4.1 Estimate of cr 

Our short distance/low T results allow an estimate for the string tension which turns out 
to be much more accurate than the existing ones at the considered value of j3. Again, as 
discussed above, we cannot use a simple minded fit to obtain a from our data. We have 
chosen to use the results of the previous analysis to extract, for each value ofR, an estimate 
of the string tension from the measurements of Q(R). The stability of the extracted values 
as a function of R also gives us a test of reliability. Based on the previous analysis of the 
corrections to the free bosonic string, we have assumed no boundary term and we have 
estimated o as 

using 63 = 0.145. The results are reported in table [HI Interestingly, we note that, for each 
one of the four values of L, the data are compatible within the errors for all the values 
R > R c , i.e. for R > 8. For each L we consider the R = 8 result as our best estimate for 
a since it is the one with the smallest uncertainty. The four values of L are independent 
experiments and we can average over them yielding the final estimate a = 0.025900(12). It 
is interesting to compare this result with the existing ones. In ref. [23J we find two different 
direct estimates obtained from two different lattice sizes: a = 0.02611(19) for L = 24 and 
a = 0.02631(13) for L = 32. Using the perturbative formula |40|) one obtains a value 
inbetween, but with a much larger uncertainty o = 0.02622(45). Our result is compatible 
with this last estimate, but it is definitely smaller than the two direct estimates. Finally 
it is also interesting to observe that the same pattern seems to be present in the results 
of [12] at (3 = 10. At this value of p, eq.flU) gives a = 0.02085(37) while ref. [12] obtains 
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o = 0.02065(6) which, as in our case, is smaller than the extrapolated one even if compatible 
within the error bars. 



R 


L = 42 


L = 48 


L = 54 


L = 60 


2 


0.031400(11) 


0.0313938(87) 


0.0313751(10) 


0.0313932(73) 


3 


0.026578(14) 


0.026569(11) 


0.026545(13) 


0.0265706(96) 


4 


0.025993(17) 


0.025983(14) 


0.025955(16) 


0.025987(12) 


5 


0.025915(20) 


0.025904(16) 


0.025873(19) 


0.025911(14) 


6 


0.025911(22) 


0.025898(18) 


0.025865(21) 


0.025909(16) 


7 


0.025914(25) 


0.025901(20) 


0.025865(24) 


0.025915(17) 


8 


0.025914(27) 


0.025902(22) 


0.025863(26) 


0.025920(19) 


9 


0.025912(29) 


0.025901(23) 


0.025859(29) 


0.025922(22) 


10 


0.025907(31) 


0.025899(25) 


0.025855(32) 


0.025922(24) 


11 


0.025902(34) 


0.025895(28) 


0.025852(35) 


0.025920(28) 


12 


0.025896(36) 


0.025893(32) 


0.025851(40) 


0.025918(33) 


13 


0.025890(40) 


0.025894(38) 


0.025848(48) 


0.025910(40) 



Table 6: Values of a extracted using eq. A52\) and the data of tabled as input. 



4.2 Scaling behaviour of 73 

In order to test the scaling behaviour of the corrections we have found at (3 = 9, we have 
performed the same analysis on a data set collected at (3 = 7.5 on a 32 2 x 48 lattice. We 
have observed the same behaviour as at (3 = 9 and - including in this case also the values 
at R = 6 and 7 since R r = 5 - we have obtained 



0.08(2) 



o 



0.03856(4) 



(53) 



From these data we estimate 73 = 0.0031(7). This value is compatible with the one ob- 
tained at j3 = 9, indicating that the observed corrections are not due to scaling violations 1 . 
This is well exemplified in the figure El where the data for A(R) for the two samples 
(/? = 9, L = 48 and (3 = 7.5, L = 48) are compared as a function of the scaling variable 
R^/a. In table [7| we report the results for Q(R) and c(R) we have obtained at j3 — 7.5. 

1 A similar analysis of possible scaling violations was performed in [21] for the Z(2) gauge theory: the 
result is consistent with what we observe here. 
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R 


Q(R) 


c(R) 


2 


0.054827(13) 


— 


3 


0.047808(17) 


0.077682(44) 


4 


0.044476(20) 


0.094185(99) 


5 


0.042658(22) 


0.10532(19) 


6 


0.041559(25) 


0.11301(33) 


7 


0.040845(27) 


0.11832(56) 


8 


0.040357(30) 


0.12201(94) 


9 


0.040008(32) 


0.1247(16) 


10 


0.039746(36) 


0.1288(28) 


11 


0.039550(41) 


0.1288(54) 


12 


0.039404(49) 


0.125(11) 


13 


0.039304(66) 


0.109(28) 



Table 7: Values of Q(R) and c(R) for the simulation at /3 = 7.5 on a 32 2 x 48 lattice. 

5 Analysis of the long distance/high T data 

In this section we discuss the results we have obtained in the long distance/high T regime. 
Our numerical simulations have been performed on a lattice with temporal extent L — 8. 
The values we have obtained for Q(R) and c(R) are reported in table [HI 

These data show that in the short distance/high T regime the behaviour of c(R) is 
completely different from the one observed in the long distance/low T case. There is no 
evidence of a convergence toward the coefficient of the Luscher term in 3 dimensions, vr/24, 
but instead a linearly rising behaviour is observed (see figure EJ. The reason is that most of 
the data are in the 2R > L regime where a linearly rising behaviour for c(R) is exactly what 
one expects from the free bosonic effective string model eq. (|3T|) . Let us stress that this 
agreement is highly non trivial and represents a very important test of the effective string 
description for the static quark potential. As a matter of fact, since eq. (13TJ) is obtained 
by modular transforming the Dedekind function, its agreement with the numerical data 
is actually a test of the whole functional form of the effective string correction 2 . In other 
words the agreement is obtained only keeping into account all the exponentially decreasing 
terms contained in the Dedekind function (which are mandatory ingredients in the modular 
transformation). Thus, even if in a rather indirect way, it could be considered as a test 
of the fact that the spectrum of the model in the short distance/low T regime should be 
the one predicted by the bosonic effective string model. This seems to disagree with the 
results recently discussed in [T3], where the authors investigate the spectrum of the stable 
excitations of the confining flux tube. From this study it turns out that below Ifm - i.e. 

2 The same holds also for the Eisenstcin functions and thus for the Nambu-Goto contribution. 
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Figure 2: Comparison of the results for A(R) as a function of the scaling variable R\fo 
for the two samples: (3 — 9, L = 48 and (3 = 7.5, L = 48 

for a range of distances where the Liischer term is already observed in the ground state 
of the static quark potential - the spectrum of the flux tube excitations is instead grossly 
distorted w.r.t. the expectations based on the effective string description. Further studies 
are needed to better understand the connection between our indirect test and the direct 
numerical evidence of |1HJ . 

Following the above comment, it is also important to notice that in the long dis- 
tance/high T regime any attempt to fit the data with a "Liischer" -type correction c(R) ~ 
7r/24 leads to a complete disagreement with the numerical simulations (see the continuous 
line in figure EJ) already for values of R comparable with the inverse temperature. This is a 
rather important remark, since this fitting procedure is very common in the literature and, 
in this regime, could lead to an apparent and erroneous rejection of the effective string 
model or to a systematic error in the estimate of o. 

The difference between the free bosonic string and the Nambu-Goto model is, in gen- 
eral, very small, but it increases as the temperature increases. The choice of L = 8 was 
motivated exactly by this consideration since at this temperature the difference between 
the two effective string predictions is definitively larger than the errors of our simulation 
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R 


Q(R) 


c(R) 


2 


0.037412(41) 


— 


3 


0.030931(49) 


0.071497(73) 


4 


0.027558(57) 


0.09488(18) 


5 


0.025497(65) 


0.11847(35) 


6 


0.024094(72) 


0.14288(60) 


7 


0.023066(80) 


0.16832(94) 


8 


0.022274(87) 


0.1948(14) 


9 


0.021643(94) 


0.2222(20) 


10 


0.02113(10) 


0.2505(27) 


11 


0.02069(11) 


0.2792(36) 


12 


0.02032(11) 


0.3080(47) 


13 


0.02000(12) 


0.3380(61) 


14 


0.01973(13) 


0.3683(76) 


15 


0.01948(14) 


0.3981(94) 


16 


0.01926(15) 


0.430(11) 


17 


0.01905(15) 


0.459(14) 


18 


0.01887(16) 


0.491(17) 


19 


0.01869(17) 


0.526(20) 



Table 8: Values of Q(R) and c(R) for L = 8 at (3 = 9. 



and we can discriminate between them. As figure 0] shows, the situation in the long dis- 
tance/high T regime is completely different from the short distance/low T one discussed 
in the previous section. In fact A(R) has now the same sign of the Nambu-Goto prediction 
and its shape and numerical value are in good agreement with it. This is one of our main 
results and suggests that, contrary to what happens in the short distance/low T regime, 
the Nambu-Goto string can give a reliable description for the static quark potential in the 
large distance/high T regime. Interestingly, a similar scenario also occurs in the 3-d Z(2) 
gauge theory [TT] . 

In principle one could interpret the deviation from the free bosonic string behaviour 
in figure 0] as due to a boundary term instead of a Nambu-Goto like correction. This 
requires b to have a value b ~ —0.8 which agrees in sign with the one that we find under 
the same assumptions in the short distance/low T regime, but it is from three to eight 
times bigger. Again, this could suggest that the deviations in A(R) might not be due to 
a boundary term. As a matter of fact the value of b could be fixed without ambiguities 
looking at different values of L since the boundary correction and the Nambu-Goto one 
have a different L dependence. However this is not easy since it would require to keep the 
high precision of our results also for larger values of L. This type of analysis was indeed 
performed in the case of the 3-d Z(2) gauge theory [21] where, thanks to the different 
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Figure 3: Comparison of the numerical results for c(R) at L — 8, /3 — 9, with the expec- 
tations for the free bosonic string (dashed-dotted line) and with the Nambu-Goto string 
model truncated at the second order (dashed line). The continuous line is what one would 
find assuming a correction to the linear behaviour like the Liischer term tt/24R. 

nature of the simulation algorithm fI7\ I2HI ITT] , this level of precision can be reached for 
arbitrary values of L and R. We plan to address this type of analysis for the 577(2) model 
in a forthcoming publication. 

Further interesting informations can be obtained from Q(R) (see figure EJ). As we 
mentioned in section ^ the behaviour of Q(R) is dominated by a. However in the present 
case we can use the very precise value for sigma obtained in section 0] In figure O we 
plot the difference (Q(R) — o) as a function of R. Thanks to the accurate estimate of 
cr, the systematic uncertainty in this difference due to a turns out to be negligible w.r.t. 
the statistical errors of Q(R). Figure El shows that the numerical estimates lie below the 
Nambu-Goto prediction truncated at the second order. On one hand, this means that 
the Nambu-Goto prediction truncated at the second order describes the data better than 
the free bosonic string, on the other hand it also suggests that the higher order terms in 
the expansion of the Nambu-Goto action could fill the remaining gap between numerical 
data and theoretical estimates. Finally let us notice that a similar agreement with the 
Nambu-Goto model was also recently observed in the large distance/high T data for the 
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Figure 4: Deviation of the results of our numerical simulation at L = 8 and of the Nambu- 
Goto string model truncated at the second order (continuous line) from the free bosonic 
string approximation. 

3-d Z(2) gauge theory. 



6 Comparison with Z(2) and SU(3) gauge theories 



It is very interesting to compare our short distance/low T results with those recently 
obtained for the Z(2) [HI HI EH and SU(3) |] gauge theories. Among the SU(3) data 
published in jH] we have selected, for our comparison, the sample at (3 = 20 (see the data 
in table 3 of ref. |Hj) which is the one with the larger set of data. For the Z(2) gauge theory 
we have chosen the sample at /3 = 0.73107 and L = 80 of ref. [21] which corresponds to 
a very low temperature T/T c = 1/20 and to a similar value for the lattice spacing. For 
the SU(2) case we have used for the comparison only the L = 42 sample which is the 
most precise one. The other three samples would add no further information since they 
are compatible within the errors. For distances R large enough, the three theories show 
the same asymptotic behaviour. Our aim is to study the higher order corrections which 
characterize the approach to the asymptotic free bosonic string limit. 
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Figure 5: Comparison of the numerical results for (Q(R) — a) at L = 8, (3 = 9, with the 
expectations for the free bosonic string (continuous line) and for the Nambu-Goto string 
truncated at the second order (dashed line). 

We plot the three data sets in figure |H1 Since the three sets correspond to different 
values of a, we have rescaled R in units of \fa in order to make a meaningful comparison. 

Let us briefly comment on figure El 

• As discussed in |Sj, the SU(3) data are rather well described by the free string 
correction only and A(R) turns out to be very small for R > R c . 

• On the contrary, both the SU(2) and the Z(2) data show rather strong deviations 
from the free string behaviour. These deviations are stronger in the Z(2) case than 
in the SU(2) one. 

• For all the three models the observed corrections definitely disagree with the Nambu- 
Goto prediction (which is not reported in the figure but it can be seen in figure [Q 
which has the opposite sign and is larger in magnitude. 

Our conclusion is that, in 3 dimensions, the terms subleading to the free bosonic string 
behaviour in Z(2), SU(2) and £77(3) gauge theories do not show any evidence of common 
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Figure 6: Values of A(R) in 3-d for the 577(3) Yang-Mills theory (diamonds), the Z(2) 
gauge theory (circles) and the sample with L = 42 of the SU(2) Yang-Mills theory 
(squares). The vertical line corresponds to R — R c — a/1.5/ct. Here A(R) is evaluated 
considering R and not R. 

features. This confirms the expectation that these subleading corrections depend on the 
string dynamics and are thus different for the different gauge groups. 

A similar comparison between the 3-d SU(2) and Z(2) gauge theories for the ground 
state energy and for the first excited string level was recently performed in [T3*l ll4|. In that 
case the two models were shown to have similar behaviours after readjusting the ratio of 
the string tension to the glueball mass in the Z(2) gauge theory. It would be interesting 
to understand if this agreement is a consequence of the rescaling or if it is instead related 
to the fact that the two models share the same center. This is one of the reasons for which 
we plan to extend our analysis also to the gauge theory with gauge group Sp(2) [2H I3U] 
which has the same center as SU(2) and Z(2). 
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7 Concluding remarks 

In this paper we have studied the corrections to the linear behaviour of the static quark 
potential in the SU(2) Yang-Mills theory in (2+1) dimensions. We have considered both 
the zero and the finite temperature regimes comparing the numerical results with the 
expectations coming from effective string descriptions of the quark interaction. Finally we 
have checked for common features among the effective string descriptions of Z(2), SU(2) 
and SU(3) gauge theories. Our results can be summarized as follows. 

• Our numerical simulations strongly support the conjecture that the ground state of 
the static quark potential in the (2+1) dimensional SU{2) Yang-Mills theory (as well 
as in all the other pure gauge theories studied up to now) is very well described by an 
effective string theory. Apart from subleading corrections - which can be observed 
at short distances and/or at high T - this effective string theory is the bosonic string 
theory originally suggested by Liischer, Symanzik and Weisz in JI]. 

• We observe deviations from the free string behaviour also in the short distance/low 
T regime. These deviations definitely disagree with the Nambu-Goto predictions. A 
similar disagreement was recently observed also in the 3-d Z(2) gauge theory [TUJ[T3J 



• In the large distance/high T regime, we clearly observe deviations from the free 
bosonic string approximation. These deviations qualitatively - and, in part, also 
quantitatively - agree with the predictions of the Nambu-Goto action truncated at 
the second perturbative order. 

• The large distance/high T regime can be obtained by a modular transformation 
of the effective string results in the short distance/low T regime. The observed 
agreement with the effective string expectation not only confirms the presence of the 
Liischer term but also suggests an indirect check of the whole functional form and, 
in particular, of the spectrum of the excited states. However we emphasize that this 
is still an open question that deserves further investigations. 

• Comparing the Z(2), SU{2) and SU(3) gauge theories, we do not find any clue for 
common features in the corrections to the Liischer term among them. 

In view of this last point, it would be important to understand if these analogies and 
differences are related to the gauge group, to its center, or are non-universal features, 
perhaps related to different values of the boundary term b. To this end it would be very 
interesting to extend this type of analysis to a wider set of pure gauge theories. For 
instance, one can consider the Yang-Mills theory with gauge group Sp(2) J23 HID] which 
has the same center Z(2) as the group SU{2) studied in this paper. 
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